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Abstract 

We consider zero-temperature, stochastic Ising models ci* with nearest-neighbor 
interactions and an initial spin configuration chosen from a symmetric Bernoulli 
distribution (corresponding physically to a deep quench). Whether a°° exists, i.e., 
whether each spin flips only finitely many times as i — > cx3 (for almost every and 
realization of the dynamics), or if not, whether every spin — or only a fraction strictly 
less than one — flips infinitely often, depends on the nature of the couplings, the 
dimension, and the lattice type. We review results, examine open questions, and 
discuss related topics. 

1 Introduction 

The behavior of different kinds of magnetic systems following a deep quench comprises a cen- 
tral topic in the study of their nonequilibrium dynamics. Rigorous and nonrigorous results 
have been obtained on different questions arising naturally in this context: the formation 
of domains, their subsequent evolution, spatial and temporal scaling properties, and related 
questions (for a review, see Ref. |1[); the persistence properties at zero and positive tempera- 
ture [0, 1^, D, ; the observed aging phenomena in both disordered and ordered systems 



(see, e.g., ^ ^ |r^); and many others. 

In this paper we will summarize results on a very basic question, whose answer is not only 
relevant to the questions mentioned above but often naturally precedes them. Put informally, 
consider a quench from infinite temperature to zero temperature, and let the system then 
evolve using standard Glauber dynamics. Will the spin configuration eventually settle down 
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to a final state, or will it continue to evolve forever (and if so, in what sense)? Here we will 
concern ourselves mostly with this question of approach to a final state in different Ising spin 
models, and will not address, except briefly in the last section, the properties of such final 
states. We now state the problem more precisely. 

Consider the stochastic process a* = cr*(a;) corresponding to the zero-temperature limit 
of Glauber dynamics for an Ising model with Hamiltonian, 

'H = - Jx,y<yx<yy ■ (1) 

{^,y} 

\\x-y\\=l 

Here 1 1 ■ 1 1 denotes Euclidean length. For now a* takes values in 5 = { — 1, +1}'^'', the space of 
(infinite-volume) spin configurations on the d- dimensional hypercubic lattice, but later (see 
Subsec. ^A\j we will examine other types of lattices as well. The initial spin configuration cr° is 
chosen from a symmetric Bernoulli product measure (denoted Pao), corresponding physically 
to a deep quench. If the spin model is disordered, then the transition rates depend on a 
realization J' of the (i.i.d., unless otherwise specified) random couplings J^^y, with (common, 
unless otherwise specified) distribution fi on R. The (continuous time) dynamics is given 
by independent (rate 1) Poisson processes at each x corresponding to those times t [think 
of these as clock rings at x] when a spin flip (o"*"'"° = — (T*~°) is considered. If the resulting 
change in energy is negative (or zero or positive), then the flip is done with probability 1 (or 
1/2 [determined, say, by a fair coin toss], or 0). We denote by the probability distribution 
on the realizations u of the dynamics. For most of this paper we consider only this single- 
spin-fiip dynamics, but in the last section will briefly discuss multi-spin-fiip dynamics. The 
joint distribution of J', 0"°, and u will be denoted P. 

A natural question in both the disordered and non- disordered models is whether cr* has 
a limit (with P-probability one) as t — >■ cx) or equivalently whether for every x, cr* flips only 
finitely many times. More generally, one may call such an x an jF-site (JF for finite) and 
otherwise an X-site (T for infinite). By translation-ergodicity, the collection of JF-sites (resp., 
X-sites) has (with P-probability one) a well-defined non-random spatial density Pjf (resp., 
Px). The densities pjr and pj depend only on d, p, and possibly lattice type, and of course 
satisfy pr + Pi = 1- We characterize the triplet {d,p, lattice type) as being type JF or X or 
Jvi (for mixed) according to whether pjr = 1 or pi = 1 or < p^r, px < 1. 

In the remainder of the paper, we review results for different different lattices, and 
a number of important special cases of p. In Sec. 0, we present results on one-dimensional 
chains for both homogeneous and disordered systems. In Sec. ^ we consider the homogeneous 
ferromagnet (or antiferromagnet) on the square lattice 7? and show it is type X. In Sec. ^ we 
review models with continuous disorder, and discuss a theorem whose consequence is that 
most such systems of interest, including ordinary random ferromagnets and the Edwards- 
Anderson (EA) spin glass |T^, are type JF for any d and lattice type. We also show why 



this theorem implies that homogeneous ferro- and antiferromagnets, on lattices (in any d) 
where each site has an odd number of neighbors, are now type J^. In Sec. |^, we discuss ±J 
spin glasses on 7? and other models with noncontinuous disorder on Z*^ that are type }A. In 



2 



Sec. ^ we summarize our findings and list a number of open problems. We also discuss there 
several situations not discussed in the bulk of the paper, including positive temperature, 
multi-spin-fiip dynamics and persistence. 



2 One-dimensional models 

In one dimension the analysis is particularly simple, and it is not hard to show that: a) 
when the couplings all have the same magnitude (regardless of sign, since in one dimension 
all such models can be gauge-transformed to the uniform ferromagnet), the model is type X; 
b) disordered models with n continuous are type JF; and c) models in which the couplings can 
(with positive probability) take on two or more discrete magnitudes are type Ai. A proof 
of the first two claims may be found in [1^, and a proof of c) in [14|. Here we summarize 



the proof of b); the proofs of a) and c) are simplifications of the corresponding d = 2 proofs 
given in Sees. |^ and H, respectively, and so will not be given separately. Modified arguments 
can be used to examine the dynamical behavior for /i's with both a continuous part and a 



single discrete magnitude; see [|1J] for details. 

A sketch of the proof of b) is as follows. Consider a site x such that | Jj^j^+il strictly 
exceeds the two neighboring coupling magnitudes; because is continuous, such sites occur 
with positive density. It is clear that once ax(Tx+i = sgn(J^^a;+i) (either initially in o"°, or 
through a subsequent flip of one of the two spins as determined by a;), neither spin will flip 
again, demonstrating already that pjr > 0. By translation-ergodicity of P, there will be 
(with P-probability one) a doubly infinite sequence of such sites Xn (with positive density). 

Consider now the interval {xn-i + 1, Xn-i + 2, . . . , By the preceding argument, there 
will be some time after which, both o"^.„_i+i and a^^ cease to flip. After that time we have a 
Markov process (restricted to the interval) with a finite state space. Because is continuous, 
each flip within the interval will strictly lower the energy, which (for fixed J) is bounded 
below, by some minimal amount. The process must therefore eventually reach an absorbing 
state in which all spins have stopped flipping. Because this argument applies to every such 
interval, a continuous /x in one dimension is type T. 



3 Homogeneous ferromagnet on 



In [13 1 it was shown that the homogeneous ferromagnet on 7? is type X; we will sketch the 
argument here. Essentially the same argument holds for the homogeneous antiferromagnet 
on 7?. In this section P refers to the joint distribution of a" and uj. 

To begin, we note that two possible absorbing states are the two uniform spin configura- 
tions: cr^ = +1 for every x and = —1 for every x. Because of global spin-flip symmetry, 
these outcomes must have equal probability p (< 1/2). But because of translation-invariance 
and translation-ergodicity of P, p must be either or 1; therefore, p = 0. The only other ab- 
sorbing states are those with one or more parallel domain walls separating strips of uniform 
+1 and —1 spin configurations; if the state is to be absorbing, then all these domain walls 
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must be parallel to either the x- or y-axis. Using the argument above, but with spin-flip 
symmetry replaced by invariance with respect to rotations by 11/2, we conclude similarly 
that each of these two sets of outcomes must also have zero probability. This shows at least 
that px > 0. 

To show that the model is type X, we need to show that every spin flips infinitely often. 
By translation- invariance and spin-flip symmetry, if pjr > 0, then the following must occur 
with positive P-probability: for some x = {xi, X2), y = {x[, X2) and z = (xi, Xg) with xi < x[ 
and X2 < X2 and for some t', a* = -|-1 and cr* = —1 and a* = —1 for all t > t' . But this 
would require that 

inf P^(a*+i i S"\a' = a) = 0, (2) 

where S" is the set of spin configurations on 7? with the values +1, —1, —1 at the sites x, ?/, z. 

But this is not so as can be seen from the following argument. Let us define w = {x'l, X2) 
and TZ to be the rectangle with corners at x,y,z and w. Any spin configuration in S" 
must have a domain wall that is contained within TZ and that either (a) connects the [x, y] 
straight line segment to the [x, z] segment or else (b) connects the [x, y] segment to the [y, w] 
segment or else (c) connects the [x, segment to the [z,w] segment. In case (a), there is 
some sequence of clock rings (and absence of rings) and coin toss outcomes within TZ during 
the time interval [t, t -|- 1] (that occurs with P^-probability bounded away from zero) that 
will move the domain wall towards the Southwest so that 0"*+^ = — 1. Similarly, in cases (b) 
or (c), the domain wall can move to the Southeast or to the Northwest so that cr*+^ = +1 or 
= +1. 



4 Models with continuous disorder 

A central result is that in any dimension (and on any lattice), models with continuous 
disorder distribution /i of finite mean are type JF. These models include EA spin glasses 
with a Gaussian p and random ferromagnets with a uniform /i. The idea behind the proof 
was already used in Sec. |^ as part of the proof that one- dimensional models with continuous 
/i are type T. Here we sketch the more general proof; for further details, see [pIS] . 
Let cr* be the value of at time t for fixed tu, and J . Let 

E{t) = -{1/2) J^y^i^l^ (3) 

y:\\x-y\\=l 

where the bar indicates an average over P. By translation-ergodicity of P, and using the 
assumption that \Jxy\ < 00, it follows that E{t) exists, is independent of x, and equals the 
energy density (i.e., the average energy per site) at time t in almost every realization of JT", 
0"°, and uj. 

Because every spin flip lowers the energy, E{t) monotonically decreases in time (note that 
-E'(O) = 0) and has a finite limit E{oc) (> —d\Jxy\). Now choose any fixed number e > 0, 
and let be the number of spin flips (over all time) of the spin at x that lower the energy 
by an amount e or greater. Then —00 < E{oo) < — eiVj so that for every x and e > 0, iVJ is 
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finite. Let be the minimum energy (magnitude) change resulting from a flip of cTx'-, then 
although varies (differently in each J) with x, it is sufficient that it is strictly positive. 

It is implicit in the proof of this result that, even without the continuity assumption on 
the distribution of couplings, with P-probability one there can be only finitely many flips 
of any spin that cause a nonzero energy change. When is continuous, the probability of 
a "tie" in any sum or difference of a given spin's nearest- neighbor coupling strengths (and 
therefore the probability of a spin flip costing zero energy) is zero, yielding the result that 
these systems are type JF. 



4.1 Homogeneous ferromagnets on lattices other than Z 

It follows from the argument above that this type result applies also to homogeneous 
Ising spin systems on lattices with an odd number of nearest neighbors, so that ties in 
energy cannot occur. (It also applies to ±J models (as defined below) on these lattices.) 
Such lattices include the two-dimensional hexagonal (or honeycomb) lattice, and the double- 
layered cubic lattices 7/^ x {0, 1} (i.e., a "ladder" when c? = 1, two horizontal planes separated 
by unit vertical distance when d = 2, and so on). 



4.2 Continuous disorder with infinite mean 

What about models where the disorder is continuous but the mean is infinite? We can show 
that a restricted class of these models are type JF; these are models where influence perco- 
lation [|I5| does not occur. Such models include all one-dimensional models with continuous 
disorder (regardless of whether the mean is finite), and strongly and highly disordered models 
0, 113. We refer the reader to [IT3 for a discussion of this situation. 



5 =bj and related models 

We now turn to models on Z"^ where 

/i = aSj^ + (1 - cy)5j2 , (4) 

with < a < 1 and Ji ^ J2. When Ji = —J2 7^ 0, we call this a ±J model. (In much of 
the literature, ±J spin glasses refer to the specific case a = 1/2.) These models have been 
analyzed in fl^, along with modifications where, e.g., /i consists of both a continuous and 



a discrete part. Because such modified models appear to be of less interest, we note here 



only that these models are generally type A4, and refer the reader to [|J] for details; for the 



remainder of this section, we confine ourselves to distributions of the form given in Eq. (H). 



The main results of |]T4| were to prove the following two assertions: a) models in which 
Ji 7^ —J2 are type Ai in any dimension, and b) ±J models in two dimensions are also type 
Ai. The first of these is much easier to prove, and we start with that case. 

We sketch here a proof for the d = 2 case; the extension to other d is straightforward. 
To show that pjr > on Z^, let |Ji| < IJ2I and consider a plaquette all of whose edges 
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have coupling J2, and all edges outside the plaquette but that connect to one of its corners 
have coupling Ji. With positive P-probability, such a coupling configuration will occur with 
the spins at the four corners initially (or eventually) all +1 (or all —1) if J2 > or else 
alternating in sign if J2 < 0. These spins will thereafter never fiip, proving that pjf > 0. 

To show that pi > 0, consider a configuration in which the couplings between some w and 
its four neighbors Zi are all Ji. Suppose also that each Zi belongs to a plaquette satisfying 
the set of conditions described in the previous paragraph, with the spins at Zi and Z2 equal 
to +1 and the spins at z^ and equal to —1. This will ensure that the spin at w fiips 
infinitely often, and because such a situation will occur with positive P-probability, px > 0. 

We now turn to a discussion of b). It is relatively easy to show that for two-dimensional 
±J models, px > 0, and we will sketch the proof of that here. It takes considerably more 
work to show that pjr > 0, and so we will present here only the idea behind the proof, and 
refer the reader to |14] for details. 

To show that pi > 0, we consider a configuration of a 5 x 5 square of sites in the dual 
lattice Z^* = + (1/2, 1/2) in which exactly 9 of the 25 sites are frustrated (corresponding 
to plaquettes in Z^ with an odd number of negative couplings): the central site Wc, two 
Southeastern sites, Wi = Wc + (1, —1) and W2 = Wc+ (2, —1); and six other sites obtained 
from {1^1,^2} by (multiple) 7r/2-rotations about Wc- Such a frustration configuration occurs 
with positive P-probability, so we are done if we can show that there always exists at least 
one Z^ site (among the 36 whose plaquettes form our 5x5 square) that has positive fiip 
rate. We now proceed to do this. 

We note first that any domain wall (i.e., a path of unsatisfied edges in Z^* [connecting a 
pair of frustrated sites]) that is not straight implies a site in Z^ with positive fiip rate; this is 
because there must exist a Z^ site with at least two unsatisfied edges. Now in our frustration 
configuration, there must be a domain wall starting from Wc- Either this domain wall already 
determines a positive fiip rate site because it is not straight, or else it runs straight out of 
the square; in the latter case, by the invariance with respect to rotations by tt/2, we may 
assume (without loss of generality) that the domain wall emanating from Wc runs to the East 
and passes just above the (dual) edge joining the two Southeastern sites, Wi and 102- But 
then there must be another domain wall starting from Wi. Either these two domain walls 
together determine a positive fiip rate site or else the second one runs from wi straight out 
of the square to the South. But then there must be a third domain wall starting from W2, 
that (together with the previous two) will determine a positive fiip rate site, no matter what 
direction it runs off to. Using as usual translation-invariance and translation-ergodicity of 
P, we conclude that pi > 0. 

A proof demonstrating that pjr > is considerably more involved, as noted, but the 
general strategy is similar. We again consider an event involving the frustration configuration 
in a finite region of Z^*, and the spin configuration in a related region of Z^. One wants to 
show that at least one of these Z^ sites will eventually have fiip rate zero and hence will fiip 
only finitely many times, thus proving p^ > 0. This is done by proving that the domain 
wall geometry in Z^* must eventually satisfy various constraints, in particular that certain 
contour events recur indefinitely with probability zero; otherwise, there would be infinitely 
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many energy-lowering flips in the fixed square with positive probability, violating the result 
presented in Sec. For further details, see fl^ . 



6 Summary and open problems 

We have studied the dynamical evolution of several categories of Ising spin models, both 
ordered and disordered, following a deep quench (from infinite to zero temperature), in all 
dimensions and for different kinds of (regular) lattices. Our concern here has centered on 
the question of existence of a final state, given the usual zero-temperature Glauber dynamics 
and the Hamiltonian (|I]). It is interesting to consider other types of situations, but these will 
be left for the future. These include, for example, other kinds of spin models (Potts, XY, 
Heisenberg) , initial spin configurations not chosen from the symmetric Bernoulli distribution, 
and others. Nevertheless, we feel that substantial progress has so far been made, and we 
review the results below, including a discussion of remaining open questions (for this type 
of model and situation). 



6.1 Review of results 

All results below are for Ising spin systems with Hamiltonian (|l]), and J ^ and lo chosen 
as discussed in Sec. |1]. A given system may have three possible dynamical outcomes: it may 
be type T, JF, or whose meanings were given in Sec. ^ 

Homogeneous ferromagnets and antiferromagnets: In one dimension and in two dimensions 
on Z^, these are type X. In any dimension on a lattice where each site has an odd number 
of nearest neighbors, they are type JF |I3| . 

Models with continuous disorder: These need to be further subdivided. The most important, 
and commonly studied, cases are those models where the coupling distribution has finite 
mean; these include ordinary (EA) spin glasses and random ferromagnets. In all dimensions 
(and for all lattices) these models are type J-" [|l3l- Another class of models that are type 
[|13| are those in which influence percolation [|1^] does not occur; these include all one- 
dimensional models with continuous disorder, and the strongly and highly disordered models 



of spin glasses and random ferromagnets |]T6|, Il7|, lIS 



±J models: In one dimension, these are type T; in two dimensions, type A4 on the square 
lattice |T^. On a lattice where each site has an odd number of neighbors, they are type 



in any dimension. 

Models with other fi: If /i is of the form a6j + PS^J + v with J>0, 0<a + /5<l, and 
V continuous and supported on [— J, +J], then it is type T in one dimension [|I^. Other 
examples are distributions of the form Eq. (|), where |Ji| 7^ IJ2I5 distributions supported 
partially on a continuous interval and partially on atoms at discrete values, and so on. These 



are type M. in all dimensions ||14|| . 
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6.2 Open problems 

For Ising spin systems with Hamiltonian (|lD, and JT, cx^ and u chosen as discussed in 
Sec. m, the problems of the dynamical outcomes of the following situations remain open: 

1) Homogeneous ferromagnets (and antiferromagnets) on the lattice Z'^ (or others except 
where each site has an odd number of neighbors) for d > 3. Numerical results 0] suggest 
that these should be type X ioi d = 3 and perhaps 4, but may possibly be type J-" {or A4) 
for d > 5. 

2) Models with continuous disorder in d > 1, where both the coupling distribution has 
infinite mean and influence percolation occurs. The proof [jl3l sketched in Sec. § already 



implies that pjr > for these models, so they are either type JF or It seems reasonable 
to conjecture that they are type JF. 

3) ± J models on Z*^ (or other lattices except where each site has an odd number of neighbors) 
for d > 3. There is little we can say about these right now, except that it would be interesting 
to relate ±J models to homogeneous ferromagnets on the same lattice. E.g., perhaps pjr(±J 
model) > Pjf (homogeneous ferromagnet). If this were the case, then demonstrating that 
homogeneous ferromagnets are type JF on Z*^ for d > 5 would immediately resolve the 
question for ±J spin glasses on those lattices. 

6.3 Related topics 

In this section we touch on several related topics. In particular, throughout this paper we 
have concerned ourselves solely with the question of convergence of ex* (with P-probability 
one) to a final state, but have not examined the properties of this final state (when it exists) 
in the different models of interest. We also have not discussed rates of convergence to the 
final state. We will briefly discuss these issues, but first will consider positive temperature. 

Positive temperatures. Here one treats the behavior of the local order parameter rather than 
that of single spins. Construction of dynamical measures, analysis of their evolution, and 
relation to pure state structure are extensively discussed in Ref. |]T^. Here we mention only 
a few relevant results. 

It should first be noted that the categorization into types X, JF, and A4 needs to be 
modified and refined at positive temperature. Without going into detail, we will simplify 
matters here by dividing systems into those where on any finite lengthscale, the system 
equilibrates (into a pure state) after a finite time (depending on (T°, u, the lengthscale, and 
J' if relevant), in the sense that interfaces cease to move across the region after that time; 
and those where this local equilibration does not occur. The latter case we call local non- 
equilibration (LNE), of which there are two types. A precise definition requires the use of a 
dynamical measure; we refer the reader to |TD| for details and to where one of the types 



[Chaotic Time Dependence) is shown to occur in a d = 1 model with disordered rates. 

A main result of [l^ is that if only a single pair, or countably many pairs (including a 
countable infinity) of pure states exists (with fixed J') , and these all have nonzero EA order 
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parameter , then LNE occurs. A corollary is that if LNE does not occur, and the limiting 
pure states have nonzero EA order parameter, then there must exist an uncountable infinity 
of pure states, with almost every pair (as the realizations of the initial state and dynamics 
vary) having overlap zero. 

One consequence of these results is that LNE (in a rough sense the positive temperature 
equivalent of type I behavior) occurs at positive temperature (with T < Tc) in the 2D 
uniform ferromagnet and (presumably) random Ising ferromagnets for d < 4. Because the 
number and structure of pure states at positive temperature in Ising spin glasses is unknown 
for d > 3 (and, from a rigorous point of view, unproved even for d = 2), occurrence of LNE 
there remains an open question. 

Multi-spin dynamics. In Ref. we examined ordinary spin glasses and random ferromag- 



nets in arbitrary dimension and considered an extension of the zero-temperature single-spin 
fiip Glauber dynamics in which rigid fiips of all lattice animals (i.e., finite connected subsets 
of Z'^, not necessarily containing the origin) up to size M spins can occur. Of particular in- 
terest here is the limit M —>■ oo. For the dynamics to remain sensible, we need to choose the 
rates for i^'-spin lattice animal fiips to decrease as K increases, so that the probability that 
any fixed spin considers a fiip in a unit time interval remains of order one, uniformly in M. 
A further requirement for the dynamics to be well-defined is that information not propagate 
arbitrarily fast throughout the lattice as M — oo. Such a dynamics can be constructed, and 
generates infinite- volume ground states 



Persistence. A topic of current interest is the P-probability p{t) at time t that a spin has 
not yet fiipped. For the homogeneous ferromagnetic Ising model on Z*^, this probability 
has been found to decay at large times as a power law p{t) ~ t~^('^) g ^ for c? < 4. 
The "persistence" exponent 6{d) is considered to be a new universal exponent governing 
nonequilibrium dynamics following a deep quench [Q. Our work J7[ shows that dependence 
on lattice type (e.g., square vs. hexagonal lattices for homogeneous ferromagnets) implies 
nonuniversality of this behavior, and moreover that the persistence phenomenon (which 
seems to require a system to be type X) is unstable to the introduction of randomness into 
the spin couplings. Moreover, in some simple systems it was shown that the decay p{t) —p{oo) 
to the final state is exponential as t ^ oo. 

Properties of the final state. We turn finally to a discussion of the properties of the limiting 
state a°° for ordinary spin glasses and random ferromagnets in any dimension. For the M- 
spin-fiip dynamics (with M finite) discussed above, the final states are energetically stable 
up to a fiip of any subset of M spins; we call these M-spin-fiip stable states. A number of 



results were obtained in [21|; a central result is that there is an uncountable infinity of these 



metastable states, in any d and for any M, and their overlap distribution is a delta-function 



at zero. For further results and discussion, we refer the reader to \21 
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